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SUMMARY 

An attempt has been made to identify possible interactions 
.between the Solar Sail and the natural plasma environment 
which might have deleterious effects on the operation of the 
sail or its associated payload. V/e have attempted to take 
conservative estimates of possible effects. In all cases, we 
have tried to report an upper bound to possible interactions. 

The net conclusion is that nothing we have foimd in the 
environment will preclude safe operation of the Solar Sail. 
However, the sail itself will perturb any natural environment 
sufficiently that observations of electric and magnetic fields 
from an a-ttached payload will not be possible. Similarly plasma 
composition and low energy particle fluxes will be severely 
affected. In some cases,, even optical observations of objects 
remote from the sail will be compromised. 

EMI induced by the plasma-sail interactions now seems to 
be much .smaller than what had originally been thought. 

Electrical-mechanical oscillations induced by the plasma 
also seem to be smaller than the original analysis had indicated. 
However, the completely rigorous treatment- of this potential 
problem is beyond the scope of this report. Ve recommend further 
study be conducted on this problem. Of particular interest 
for operation in the near earth environment will be a rigorous 
analysis of the interaction of the sail with natural plasma 
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1. INTRODUCTION 

The purpose of this report is to assess the, interaction 
"between the Solar Sail and the plasma environments which will 
"be encountered in various missions in order to identify any 
deleterious affects upon the performance of "the Solar Sail 
System. In addition we examine deployment problems due to 
static electricity since these fall naturally into a general 
analysis of the electrostatic properties of the Solar Sail. 

Our approach is both qualitative and quantitative, with an 
emphasis upon worst case engineering analysis. We are looking 
for the "show stoppers", those critical problems which can 
cause a major system failure. A wide range of problems is 
covered, in as much depth as time has allowed. Areas needing 
further study are indicated, as is our present state of 
understanding. 

We begin by examining the solar wind environment. Satellite 
data is used to develop a model of the solar wind which can then 
be used to determine the voltage to which a Solar Sail might be 

f 

expected to charge. The environmental definition is done from 
the perspective of spacecraft charging. Thus when parameter 
choices involving Judgement must be made, we shall always 
make choices between observed parameters which tend to give the 
largest value for the electron temperature, as this quantity 
determines the voltage to which the spacecraft will charge, and 
ve prefer to overestimate rather than underestimate this voltage. 
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Having defined the solar wind environment we then proceed 
to perform an analysis of the electrostatic charging of the 
Solar Sail by the solar wind. The first and most important 
problem addressed is differential charging. We examine the 
possibility of developing large potential drops across the 
various regions of the spacecraft body, with the possibility 
of high voltage breakdown. A second consideration is the 
disturbance of the plasma environment about the sail, and 
how this disturbance affects measurements of this environment. 

We next examine the use of the Solar Sail in the near 
earth environment. Again a general charging analysis is 
performed. Electromechanical effects (resonant interaction 
with PC4 oscillations) specific to this environment are 
examined. 

Electrostatic forces acting upon the sail arise from the 
gnvironment and from frictional effects or contact potentials. 
These latter are important during deployment, which represents 
a special case for analysis. An assessment of the minimum forces 
anticipated is performed, as part of a worst case analysis of 
the deployment problem. 

Another problem arises from the forces associated with 
electrostatic charging by the environment. Electromechanical 
oscillations arise from the coupling of waves on the sail to 
the surrounding plasma. We have analyzed this problem in an 
attempt to determine if the interaction with the solar wind 
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vill cause unstable oscillations of the sail. Several cases 
are examined and dimensionless parameters characterizing 
potential instabilities are developed. 

The existence of plasma flow across a magnetic field in 
the solar wind gives rise to still another force, the induction 
force. This force and a corresponding induction current are 
analyzed to determine their impact upon sail performance. The 
modification of the photo electron sheath by these forces is 
examined. 

For the Solar Sail to be useful it must transport a 
payload. We examine the effects of Solar Sail operation 
upon the payload. Areas of concern are; 1) safety, 2) effects 
on particle measurements, 3) effects on field measurements and 
4) electromagnetic interference (EMI). Early concerns with 
EMI have lessened substantially as the magnitude of the sail 
charging problem became better understood but we mention it 

for the sake of completeness. 

The need for experimental work is indicated in a final 
section. A few simple experiments might be especially useful 
in assessing the problems associated with static electrification 
during deployment. 
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2. THE SOLAR WIND ENVIRONMENT 

liiie solar wind is a supersonic predominantly electron- 
proton plasma flowing approximately radially outward from the 
sun. It is generated by the expansion of the solar corona. 
Between r — «3 snd past 1 .0 a.u. the density n varies 

approximately as 1/r2, the flow velocity 1 is approximately 
constant. These statements however describe only the time 
average behavior. In fact very large fluctuations are normal. 
They are due in part to the fact that there is continuous 
activity near the base of the corona and that the corona has 
important structure which rotates with the sun (e.g. coronal 
holes). Thus present understanding of the solar wind recognizes 
the existence of high speed particle streams emanating from 
coronal holes, and shock fronts which build up when fast plasma 
overtakes cold plasma. For this reason our environmental 
specification will rely on satellite data. The vast bulk of 
the observations are at 1.0 a.u. ^^*^’^^and there do exist 
theoretical models which allow extrapolation to .3 a.u. 

We shall examine this scaling as well as make use of the limited 

(5) 

data at .3 a.u. that we have been able to locate on Helios 
satellite observations. The quantities we require are the 
flow speed, density and electron and proton temperatures. 

Since the distribution functions are not isotropic Maxwellians 
some interpretation of experimental data is required. In the 
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following we shall always make choices between observed 
parameters that tend to give the largest value for the electron 
temperature, as this quantity determines the voltage to which 
a spacecraft will charge, and we prefer to overestimate .rather 
than underestimate this voltage. 

a. Distribution Functions 

Figure 2ol illustrates recent data on proton temperature, 
flow speed and density from 1.0 a.u. to .3 a.u. Flow velocities 
as high as 800 km/sec are evident. The average density behaves 
as n=1/r^. The peak density near .3 a.u. is about 300/cm^, 
and we shall use this as the most extreme condition for 
charging calculations 5 100/cm^ would be a good average value. 

From the proton temperature we wish to deduce a value of the 
electron temperature. The scaling calculations of Ref. 5, see 
Figure 2.2, seem to indicate about a 20 - ^0% change in electron 
temperature between 1.0 a.u. and .3 a.u. Figure 2.1 indicates a 
roughly similar change in the maximum proton temperature. We 
have not been able to locate electron temperature data from Helios 
and are not sure that it has been published yet. We shall 'assume 
that the more thoroughly investigated properties of the proton 
and electron distributions at 1.0 a.u. hold at .3 j make 
use of the most extreme electron temperatxire that can be deduced 
in this manner. Figure 2.3 shows a plot of electron and proton 
temperatures at 1.0 a.u. observed with the Vela 4 satellite. 
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Figure 2.1, Solar wind bulk velocity (heavy line), density 
(thin line), and temperature (dotted line) as 
measured by Helios A (Ref. 5). 
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Figure 2.2. 


Temperature and density variation with radius 
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ELECTRON - PROTOK TEMPERATURE , D£ NSITY . & BULK SPEED IN THE SOLAR WIND 



SOLAR ECLIPTIC COORDINATES: 
tONCITUDt *' LATITUDE -8»* 


Figure 2.3. Solar wind parameters at 1.0 a.u. (Ref. 1). 
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The maximiom and minimiim values of the proton temperature , 

T and T . , reflect the anisotropy of the distribution with 
respect to the magnetic ‘ field direction. For our purposes we 
can take the electron distribution as isotropic, and according 
to Figure 2.3 we can take the electron temperature to be the 
maximum proton temperature. The proton distribution is fairly 
well represented as an anisotropic Maxwellian. The proton 

anisotropy ratio is given by K = ^jnax'^^min a*u. 

it ranges from 1.1 to in "the data of Ref. 2. The average 
temperature that is typically constructed is defined as 

T = 1/3 = 1/3 W '1 t>- 

Figure 2.3 indicates values ’l.^xIO^V, and T^^^"5x10‘* °K, 

or T'S.BxlO*^.' This is somewhat lower than the value of "3x10^°K 

X 

apparent in Figure 2.1. We shall proceed on the assumption that 
3x10^ °K represents the average temperature T at 1.0 a.u. , 
noting that the bulk speeds observed in the Helios data are con- 
siderably higher than that shown in the Vela data, and that in- 

C2) 

creasing V is correlated with increasing proton temperature • 
Identifying T with! and using the maximum value of K we find 

’’e" ■^/K - • 

The average proton anisotropy increases with radius because the 
magnetic moment W_l/B is constant and B is decreasing. This 
being the case it appears that taking Tg~2T should be a safe 
upper limit on the electron temperature. From Figure 2.1 
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at .3 a.-u. we find then T ‘6x10^ °K at maximum, and thus 

T <1.2x10^ We thus have deduced so far the extreme 

e— 

plasma parameters 

V'\'300 - 800 km/sec 

n'^l 00-300/ cm^ 

T <1 .2x10^ °K 
e— 

T'6x105 

^max"''-9T = 1.14 x10^ °K 
W " 3.3x105 Ok. 

Several coDiments are in order. 

1 ) The proton temperature is not too important for charging 
calculations because the proton directed velocity is large 
compared to thermal velocity. 

2) The observations illustrated in Figure 2.1, 2.2 were 
made during periods of low average solar activity. We have 

not located corresponding data on active periods for comparison. 
IP any problems were found with solar wind spacecraft charging 
it would be important to make a more thorough study. 

3) Ko data is available for high solar latitude (beyond~10°) . 
Various models indicate variations by factors of order 2 

for quiet conditions (6). 


12 



The Electron Distribution 


It is important for charging calculations to recognize 
that although the electron distribution is approximately 
isotropic it is not a single Maxwellian. It can be adequately 
represented as the sum of a cold or "core" Maxwellian which 
contains the bulk of the electrons, and a hot or "halo" 
Maxwellian. Reference 3 gives values for the number of hot 
electrons to total ranging from .O38+.OI7 to 

.071+. 028 j and the temperature ratio 


Rm - 


_ T. 


H/T^ = 


8x10 5 °K 

1 .25x^0"^ 


■'K 


6.4 


We assume this ratio holds at .3 a.u. and identify with the 


average temperature 

Te = = 1.54T^ _ 

using Np/Np+N^“. 1 , the maximum value. 

We find 

T =■ T^/1.54 = 7 . 79 x 1 Q 5 °K = 67.2 eV 
Tp = = 5 x 1 0 ^ °K = 431 eV 

We shall use these to characterize the electron distribution. 

In Ref. 5 a feattire of the solar wind electrons was 
suggested. This is that the halo electrons are actually the 
scattered portion of an electron beam component of the solar 
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wind. The halo electrons appear to have roughly the same 
energy as the observed beam component. We shall not attempt 
to incorporate this new feature into o\ir charging estimates 
but work with the more conventional two temperature Maxwellian. 
Thus for the electron distribution we take 

where fjj^CT) is the normalized Maxwellian characterized by 
temperature T. In f^^ we neglect the flow speed because it is 
small compared to thermal speed for electrons. 
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3. SAIL CHARGING IN THE SOLAR WIND 

The interaction of a spacecraft with the ambient plasma 
is an important consideration for a number of reasons. The 
first and most important of these is the possibility of 
developing large potential drops across various regions -of 
the spacecraft bo.dy, with the possibility of high voltage 
breakdown. This could lead to electromagnetic interference 
with consequent electronic malfunction and possibly permanent 
damage. In the case of the solar sail breakdown occuring m 
the kapton sail film could lead to premature degradation of 
the. sail itself. A second common consideration is the disturbance 
in the ambient electric and magnetic fields and particle 
environments produced by the spacecraft, and how measurements 
of these environments are affected. To give any final answer 
to this question would require more detailed analysis than 
could be provided in the term of this study; we have provided 
a short section suggesting a few possible effects of the sail 
on payload operation. It is true that plans call for jetissoning 
the sail before encounter with Halley's comet, in which 
case the problem becomes a more conventional one. On the other 
hand the cranking orbit offers a unique opportunity to 
examine the solar wind at all latitudes; data at present exists 
for only a few degrees out of the eliptic. 

Summ^. The charging characteristics of the sail 
depend strongly on vhether or not the front (sunny) side and 
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rear (shady) side are electrically shorted together. If 
disconnected the front side should fl.oat. at some 1-10- volts 
positive with respect to space and the rear side could float 
to 300 volts negative. The shorted sail should float at 
1-10 volts positive, a safe and desirable situation. The 
small positive potentials result because of large . fluxes of 
photoelectrons that are emitted and the proton flux from the 
directed flow of solar wind. For heliogyro blade surfaces 
near parallel to the solar wind these factors are greatly 
reduced and large negative potentials (-200 to -hOO volt) 
could develop over the whole blade. We have analyzed this 
possibility and find it unlikely that this could occur. If 
all blades were shorted together it would be virtually 
impossible. 

Spacecraft charging and wake . When imbedded in a plasma 
a spacecraft floats to a potential such that the various 

currents to the spacecraft surface, which are all functions 
of if, sum to zero. The major current contributors are plasma 
electrons minus associated secondaries, plasma protons plus 
associated secondary electrons, and photoelectrons emitted 
from sunlit sTjrfaces. 

The sail is imbedded in the flowing solar wind plasma. 
The dynamics of this plasma is dominated by the flow velocity 
of the protons (which are the dominant ion constituent). The 
protons are supersonic, that is the flow velocity 
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7>>V , the proton thermal velocity (V/Vp = M'V'4-10). On 

the other hand typically v^ however we shall consider 

the possibility that in some circumstances v^ does not exceed 
V by a huge amount. Under these conditions (M>>1) a cavity 
is formed in the plasma behind the spacecraft. There are two 
major factors affecting the dynamics of this cavity. If all 
ions had precisely velocity ^ and the ones hitting the sail 
were perfectly absorbed a cylindrical cavity whose cross 
section was that of the sail frontal cross section would form, 
Figure 3.1. The ion thermal spread gives the average ion a 
velocity component v in the transverse direction. This alone 
would cause the cavity to fill in a distance downstream 
b'v(Vp/V)2, where w is the smallest transverse sail dimension. 
Aside from ion thermal effects the wake will fill in due to 
plasma electric fields. These fields come about f or t wo 
reasons, 1) the electrons being much faster than the ions 
tend to drag them into the cavity (ambipolar field), 2) the 
spacecraft surface is typically charged to some potential., which 
also influences the trajectories. 


The problems of theoretically calculating the structure 
of the disturbed plasma (frequently referred to as the wake) 
moving around a body in space requires solving a complicated 
system of coupled nonlinear partial integro-differential 
equations. The equations consist of the Vlasov (collisionless 
Boltzmann) equation for the ions and electrons, and the Poisson 
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eqiiation relating the electric field to the particle distributions. 
The final solution must give zero electric current to the 
spacecraft, which means that its surface potential must be self- 
consistently included. Since such a calculation was well beyond 
the scope of this study we have resorted to an examination of 
the literature to determine if some useful physical information 
could be extracted. We have concentrated on those factors 
which would be important for the charging of the sail. That 
is we have attempted to extract information about electric 
ciirrents to the sail independent of detailed wake morphology. 

There have been numerous calculations of the wake of a 
large spacecraft. These are reviewed in References 1-6. In 
none of these have we found any discussion of plasma properties 
exactly at the rear surface of the spacecraft. Many of the 
calculations have in fact been carried out for planar surfaces 
(in the limit of large Mach number M the calculations for an 
arbitrarily shaped body reduce to those for a disc with the 
shape of the body’ frontal cross-section).. Typically 
approximations are made -which are simply wrong near the rear 
surface. We feel it is still possible to obtain a useful 
estimate of the plasma current to the rear surface. 

The largest particle flux in the plasma is the electron 
flux, because they are travelling much faster than the protons. 

In the absence of photoemission this causes the spacecraft 
surface to come to a potential which is negative. In this 
circumstance the electron current density to a surface element 
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under quite general conditions is the sum of 



for each Maxwellian component of the electron distribution 
at infinity. This expression will be reduced on concave regions 
of the body surface since some current producing trajectories 
coming in from infinity are cut off by the body. For the solar 
sail geometry, and in particular the heliogyro, this should 
be a negligible factor. If there exists a potential well 
(electron surplus) astride the trajectories leading to the 
surface element in question this too will reduce the above 
expression since the electron potential energy function will 
have a hill. For a disc shaped spacecraft whose radius is 
less than five Debye lengths Xjj such wells apparently do not 
exist (4), although they certainly could behind the square 
sail since the square sail is some 80-250 Debye lengths wide. 

We shall be focusing on the heliogyro. For it the blade width w 
is anywhere from to depending on solar distahce and 

conditions. The work that we have examined for infinitely long 
blade like configurations (6) indicates that shallow potential 
wells may exist. These calculations may not be definitive however 
since the ion thermal velocity spread was neglected, and they 
seem hard to reconcile with the disc results of Ref. 4. Since 
potential wells reduce the electron flux^^^ we feel it is 
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prudent to use the electron current calculated as if they did 

not exist, that is This will give a worst case value 

for the negative potential on the rear of the unshorted sail. 

If the sail is positive the electron flux will be increased 

over simply the surface integral of For example if a 

positively charged plasma pro be is cy lindrical in shape its 

I ^ 


effective collection area isjl + 




kT. 


times its actual area. 


and we might expect a similar effect to hold for a bladelike 
probe. Since the positive potentials we shall find are typically 
just a few volts ef /kT^<< 1 we shall take the enhancement 
factor to be unity. This has the' effect of overstating 
small slightly positive potentials which are not serious for 
sail operations anyway. 


It is reasonable in our estimation to take the ion flux 
to the rear of the sail to be simply zero. For an uncharged 
sail it requires protons moving towards the sun to hit the rear 
of the sail, and these are a very small fraction for solar wind 
Mach numbers M. Figure 3.2 shows surfaces of constant ion 
concentration^ behind a semi-plane. The ion density is clearly 
falling rapidly as the rear surface is approached. For a 
negatively charged rear s'urface additional ions will be attracted 
by the electric field. These will tend to reduce the field. 

Thus it should be a good first approximation to neglect these 
ions. We shall use the leakage current through the kapton to 
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balance the electron current. In one sense this is academic. 

The presence of large negative rearside potentials is undesirable 
enough, and its remedy, the shorting together of front and 
rear, simple enough that it almost certainly will be done. On 
the other, hand the rear side is vulnerable to any proton fluxes 

O 

that may exist because its metallic coating is only 50-100 A 
thick. Calculations have suggested that the plasma behind 
the body is unstable.^ ^ and this could cause particles to 

be heated and deflected to the rear surface. In Figure 3.3 
are illustrated the proton distribution function as a function 
of angle from the sail edge, showing that it becomes beamlike, 
and the calculated region of instability. We do not know 
if such instability will be important. 

In summary we shall be using r for positive surfaces , 
r ) for negative surfaces, and the proton flux to be 
calculated in the following for the front of the sail. We 
shall be neglecting the potential gradient that exists along 
the sail surface due to the electric field induced by the 
VxB of the plasma (Section 7). 

Plasma Currents to Sail . The external currents we shall 
take into account are plasma protons, plasma electrons, photo- 
electrons, and a rough estimate of secondary electrons due to 
electron impact. 
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T = M tan e 
= V_ v/2 -X 
■ '^p " 



Figure 3»3(a). The ion distribution function g(u,T) for 

various values of velocity angle with respect 
to solar wind direction (Ref. 3, 7). See 
Figure 3.1 for coordinates. 


I 

p 

t =z/(w/2) 



Figure 3.3(b). The boundary of the region of instability for 

flow past a plate (Ref. .3, 7). Ion acoustic 
waves are unstable in the region above the 
curve . 
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1) Plasma Protons: The proton distribution function is 

taken as a Maxwellian with different temperatures parallel and 
perpendicular to the magnetic field t which are on the order 
of 10^ but can extend higher (see Section 2 )• We assume 
that B and the directed flow velocity lie perpendicular to 
the plane of the sail (it will turn out that this does not 
matter). The distribution function is thus 


/m \ f'^'D 


2TTkT|, 




2kT 


2kT 


Jl 


the flux to the sail is 


-I 


/m. 


dv^dvydv^.v^fp = n|^^ 


n 


■)*r 






min 


V is the minimum positive velocity a proton can have and 
min 

still reach the surface: = 0 for a surface at zero or 

negative potential and v . j, = ^ positively charged 

^T) • • 

surface. We normalize to the plasma temperature, defining 
4)=e^ /kT|j . The flux result is 

= n V iM(1 + erf(M- /7))j . 

P PL 


Regarding erf as an odd function includes the case $> M. Solar 
wind proton Mach numbers are of the order M = For zero 

body potential ^the flux is very well approximated by = nV. 

Thus the amsotropy with respect to B is insignificant. In 
Figure 3.4 we have plotted the flux normalized to as a 
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0 5000 10000 

Figure 3.4. Proton flux as a function of positive sail potential. 


function of y . VJhat this plot shows is that very high 
positive potentials are required to significantly retard the 
protons, and is a good approximation for potentials up to 
^-1000 volts. Since such high positive potentials are unlikely 
to be encountered the proton flux w,e shall use in calculations 
will be Tp = Tpp. 

2) Plasma Electrons: The integral calculated above for the 

■protons applies to the electrons. For these M<<1 typically. 

e^ /kT 

For K = 0 the electrons fluxes are and r__ e ^ , 

for each Maxwellian component. 

3) Photoelectrons: The charging of a body in space is 

strongly influenced by photoemission of electrons from its 
sunlit s'urfaces. The emitted fl'ux depends on the solar photon 
energy spectrum and the emitting material. If the body potential 

is positive with respect to that at large distance the fl'ux 
will also depend on ^ because low energy electrons will be 
attracted back to the surface. Extensive work has been carried 
out'^^^ measuring the photoelectron spectriun as a function of 
incident photon energy. These are then folded in with the solar 
spectr'um and integrated o'ver angles of emission to give the 
sunlight emission spectrum and current flow from a surface. 

/ D \ 

The total flux from aluminum at 1.0 a.u. is given as'‘ 

I = 2.6x10^ °/cm2-sec, and would be I = 2.9x10^ ^/cm^-sec at 
• 3 a.,u. The. perpendicular energy distribution, shown in 
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Figure 2.1-d, Ref. 8, can be reasonably well represented as 
the .sum .of two- -Ma-xwe-lli-ans- 

f(E) eC ( j 

with ai = 1.1, 02 = El = .625'eV, E2 = 2.0 eV. 

The total flux from a positive surface is 

■■Bh r. 


'min 


normalized to 1^, where 


in y 


2ey 


m. 


£ when y >0 and - 0 


when y <0. The final expression we find is 

3 


^ph 

for y g>0 


/ 1.86x1010(1.0 

I 2.06x1011( .3 


a.u. ) \ . [g-1.6e5j:ev) 
a.u. ) j L 





If) Secondary Electrons from Electron Impact: Secondary emission 

processes seem to have a threshold for occurence at about 
40-50 eV incident electron energy^^\ The secondary yields 
rise rapidly and maximize for Al and AI2O3 near 300 eV^ and 
at similar energies for many other materials. We have found 
no specific data for Or and its oxide,, and so shall use Al 
numbers for estimates. This means that in the quiet solar 
wind (T <v20eV) secondary emission is negligible. Under extreme 
conditions T 'wlOOeV. and secondary emission could influence the 
equilibrium potential. We shall present results both with and 
without inclusion of a secondary emission estimate. 
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A commonly used formula that represents reasonably well 
the secondary yield S(E) = number of secondaries/number of 
primaries ^ 

E'E) = 1 -^ 

Here E is the incident electron energy, is the energy of 

maximum secondary emission (»>'300eV). For A1 

S %2.6 for Al?0:i. For metallic oxides the yields can be 
max 

much higher than for the pure metal (S(E) can go as high as 10). 
For electrons with E<?00eV the penetration depth is only some 
10 atomic layers, so that the intrinsic oxide layers on Or 
and A1 might provide the higher yield factor. However the 
‘oxide layers decompose due to electron impact. Thus it is 
probably prudent to use the smaller pure metal yield to 
’calculate lover limits for negative potentials. 


Integration of S(E) over the electron distribution which 
we take to be a single Maxwellian with T^'^loOeV yields the 
secondary current^^^^ 

~ Sjnax 

This expression is convenient because it says that the c\irrent 
of secondaries is proportional to the primary current. The 
quantity in brackets is 

i^erfc(z) = ^ j ® ^ 

Using A1 metal numbers we find . 
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Although this factor is large it does not alter greatly the 
expected negative potential on the rear of an unshorted sail. 

This is because the correction enters logarithmically in the 
current balance equation. 

Unshorted Sail - Shadv Side . The proton density on the 
rear side of the sail is probably very small as discussed above. 

This being the case we shall consider other currents to balance 
with the electrons. The negative potential yielded by these 
computations could then be used to estimate the proton flux, 
and a new equilibrium computed. Neglecting the protons will 
give a worst case negative potential. On the shady side there 
is no photoemission and one must look to secondary processes 
to provide the currents that balance the electrons. Since 
the secondary current we are- using is just a numerical factor 
times the primary current it cannot be used to balance the 
primary current. Equilibria between electrons and their secondaries 
do exist if more detailed expressions for the secondaries are 
used^^^\ We shall instead consider the leakage current through 
the kapton sail film. The kapton conductivity increases with 
electric field, and this turns out to be an important process 
limiting the potential on the rear side. 

To include the effect of leakage current we model the 
sail as a parallel plate capacitor with the front plate held 
at zero potential (this is close to the actual value expected). 

We determine the potential of the rear side by balancing the 
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leakage carrent with the net electron current. According to 
one model (Ref. 9, P* 3*29) the conductivity is given by 

0 = [ 2 + cosh j mho/meter. 


E is the field in volts/meter, i.e. E = ifg(volt) /2. 54x10 meter 

Here o^CT) '^1 0” ^ ^-1 0“ ^ ^mho/meter at T'^300°K 500 K 

implies we will find an upper bound oniyg|); and 6p/2kT'v7. 8x10 
A good approximation is 

0 = 1=„(T) , 


sec” 


or in c.g.s. units 

O = (1 .5x10” 7-1 .5x10’ 
with in volts. The current density is then 

j = oE = .( 1 .97x10”^-1 .97x10"5)':f esu/cnf -sec. 


env. 


kT 


This is set equal to the electron current E fg = 2 ^=^ e 

both with and without the correction factor for secondaries. 
The resultant potentials for a variety of conditions are 


displayed in Table 3,1. 

These potentials are not particularly high compared to , 
magnetospheric charging under substorm-eclipse conditions. 

The breakdown potential of the film is of the order of 1 kV, 
which may be sufficiently above these values. On the other 
hand under extreme conditions the potentials might go higher. 
Perhaps discharge could occur for reasons other than dielectric 
breakdown, such as being induced by material damage due to 
micrometeoroids. With these voltages a large amount of energy 
could be dissipated all at once in one place. The energy 
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stored in one heliogyro blade is some 5x10**5oules. Since 
shorting the front and back side together seems to present 
no special problems this is probably advisable. 

Unshorted Sail - Snnnv Side . The major currents to the 
front side are solar vind electrons minus secondaries, photo-, 
electrons, and solar wind protons plus associated secondaries. 

It should be noted that the heliogyro will be operated in modes 
such that the blades will be at angles other than 90° to the 
spacecraft-sun line. The plasma electron current to the blade 
will not be a function of angle because the electrons are 
nearly isotropic (aside from some possible shadowing effects 
which are expected to be small because the gyroradius is 
large compared to the blade width). We shall present below, 
estimates of the potential as a function of angle in the 
context of the shorted sail. 

Using the flux expressions presented above we have 
computed the 'equilibrium potentials displayed in Table 3.2. 

The potentials come out positive so that is appropriate.. 

At .3 a.u. we have used the highest electron temperature 
Tg'v-IOOeV (Tp = 6. lev, Tj^ = 431 eV). The results both with and 
without the estimate of electron secondaries are given. We 
note that in some cases with secondaries included the corrected 
electron flux is too small to balance the directed pr.oton flux. 
We have indicated these cases with an asterisk (*), and have 
included in parenthesis revised flux and potential estimates 
that take into account the directed electron velocity component. 
In some cases this correction is sufficient to yield an 
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1.0 a.u. Quiet Condit ion 

n / 2kT< 
= 20eV = p J=- I : 


n(cm“3) V(km, 

10 3x1 0 ’ 

10 8x10’ 

a) No Secondaries 


7-5x10 


1 .22x10 


8x10 


-3 a.u. 


T = lOOeV 


100^ 300 2.8x10^3-1x109) 3x10^ * (-H3) 

100 800 2.8x109(3.3x109) 8x1 p9 * ( * ) 

Too 300 8.4x109(9.24x109) 9 x1Q9 |* (+ 12) 

300 800 8.4x109(9.9x109; 2.4x10i 


ii 


b) Secondaries included 


.3 a.u. 

T. = lOOeT 

w 


= 67eV 


Tg = 431 eV 


1 .6x1019 


1.6x1019 


3x109 


8x109 


300 300 4.75x1019 9x1 Q9 


300 I 800 4.75x10^9 2.4 x1Q19 |+2.75 

c) No Secondaries 

Table 3.2, Sunny side potentials of unshorted sail 
Normal sunlight and proton incidence. 
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eqmlibri-um. "What will happen is that the secondaries which 
have only a few volts energy will be cut off as the potential 
rises. The potential cannot rise above that which will cut 
off the secondaries because the results with no secondaries 
included indicate that equilibrium is attained. It does seem 
remotely possible however that if the proton density were 
substantially higher than the electron density for a period 
of time (see for example Figure 3,3.) the proton flux could 
exceed the electron flux and give a fairly large positive 
voltage; the decrease of the proton flux with voltage and 
the increase of the effective area for electron collection 
are slowly varying functions of voltage. 

Shorted Sail - Normal Incidence . For the shorted sail- 
we still find all positive voltage equilibria. This means 
that the total plasma electron current is computed by 
doubling the entries for the unshorted sunny side. The 
results are presented in Table 3,3, Again a few cases are 
found where the electrons are overbalanced, and the above 
considerations apply. 

Shorted Sail - Oblique Incidence . The heliogyro will 

be operated in modes where a nimber of blades could lie near 

parallel to the spacecraft-sun line, Figure 3.5, In the limit 

0=0, assuming a perfectly flat blade, the photoemission 

goes to zero and the proton flux becomes the thermal rather 

than the directed flux. In this position' large negative 

potentials can arise. In the following we construct a 

simplified model of the current equilibrium as a function of 
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1.0 a.u. Quiet Conditions 


Tg = 20eV Tg - 2 /TtJ-^ 

_ _ € _ 


V(km/sec) r^Ccm'^sec" M 


1.5x10® 3x10' 


1.5x10® 8x10' 


a) No Secondaries included 


(volts) 


+ 3.75 


+ 4,7 


■3 a.u. 


n(.177y. /2kT_ . 

re 


T = lOOeV, T„ = 67eT T„ = 431 eV 



100 poo 


100 800 


5.6x10 9 


5.6x10®(6.6x10®) 8x10® * (*) 


1.68x10^® 




xIO® 


300 300 

300 800 1 .68x10i®(2.11x10i®)| 2.4x10^ 


b) Secondaries Included 


• 3 1 

T := lOOeV = 67eV T„ = 431 eY 

e 0 — _ 


II 


100 300 3.16x10^® 


100 800 3.16x10^° 


9-5x101® 


300 30 


300 800 


3x10® 


8x10® 


9x10® 


2.4x101® 


+2.45 


+ 2.7 


+1.05 


c) Secondaries not included 

Table 3.3. Shorted Sail Potentials. Nornsal sunlight 
and proton incidence 





































the blade angle e. The conclusion is that it is very unlikely 
that the heliogyro could ever attain large negative potentials, 
due to angular positioning. Photoemission is always sufficient 
to hold the sail a few volts positive. 


To construct a simple model we neglect thermal motion 
along the flow velocity direction (z) and along the blade length 
(y). We also neglect the fluxes to the rear of the blade. 

This will not significantly affect the result. The proton 
flux is 

r = f (v) v'hd^v, 

P P 

where 

h = 2 cos e' -X sin e' , 

V = v^x + v^z, 

v*n = v^, cosQ'-v^ sin e’, 

= “(lie)* e- af V 5(v^-V). 


Inserting these expressions we find 

"V tane 


rp(6) 




dv (V sin e -v„ cose )e 

X 


a— V 2 
2kT ^x 


or 


[ 1 n / 2kT^ 

1+erf(M tan e)j + 2 — v 

* P 


-(M tan e)^ 


It is easy to check that this formula has the correct limiting 
behavior in the limits e-^0, 90°. To model the photoemission 

we simply use the previously calculated flux multiplied by 
sin 6 . For an example we have taken n = 100, photoemission 
at .3 a.u., and T^ = Tp = 80eV. We find ifg(e=0)==-300 volt, 
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^ (e = tt/2) = +5 volt. 3g(e) is plotted in Figure 3,6. We 
see that photoemission holds the blade positive unless e is 
very close to zero. The angle e varies significantly along 
the blade for the real heliogyro. Thus it is unlikely that 
lar.ge negative potentials will ever occur. The reason, for 
this behavior of y_(e) is the large value of the primary 
photoemission flux. At 6= 90 °, iPg = +5 volt, some 98 ^ of 
the photoelectrons are attracted back to the sail. Thus the 
primary flux can be very small (sin 6<<1) and Sg becomes 
less positive. If the blade were built in electrically isolated 
segments a voltage differential might arise. 


• Proton Flux to Rear of Sail . The above calculation for 
j’g(e) can be adapted to estimate the proton flux to the rear 
of the sail. This is accomplished by letting fi— »-h and 


integrating over the range V tan e< v <“. We find the expression 

sine(1-erf(M tane))+ cosee’^^ tane) j ^ 


_ n 

" 2 IT 




rj^(90°) = 0 . 


The ratio rT,(6)/r is plotted in Figure 3.7, We see that unless 
the sail is lying almost parallel to the spacecraft-sun line 
is a very small fraction of (.0065 at ^-5°); fj, itself 

is only a fraction 1/2 JUli = ,043 (M=6.5) the directed flux nV. 
At 45° is a fraction 3x1 0’*^ of the directed flux. It 
vould appear that even though the metallic coating on the rear 
is only some 5-1 OJS as thick as the VDA on the front the 

proton deposition in the kapton represented by the rearslde 
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fliix should be negligible compared to that coming in from 
the front. 
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4. USE OF SOLAR SAIL IN NEAR EARTH ENVIRONMENT 


Several possible missions have been proposed for 
the sail in the vicinity of earth. Examples are (1) 

a) Halley's Comet and Mars Sample and Beturn missions 
deployment at low altitude (“800-1000 km) 
followed by solar radiation pressure driven 
orbit raising. 

b) Geosynchronous orbit (GSO) operations; Use of' 
solar sail for stationkeeping in and orbital 
transfer along GSO. 

c) Atmospheric braking of returning interplanetary 
mission. 

We shall discuss in the following some possible problems 
with earth orbit sail operations in the area of spacecraft 
charging and electromechanical effects within the magneto- 
sphere. No problems which would prevent earth orbit 
operations have been uncovered. 

A major constraint on near earth missions is the low 
thrust to mass ratio characteristic of the solar sail. For 
example in the Halley's Comet configuration only about 20% 
of the mass is in the payload. One finds that it would 
take on the order of six months to raise the orbit from 
1000 km to GSO and perhaps another two months for release 
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from earth orbit (800 km is about the minimxm altitude 
a-t which the sail could be operated, the atmospheric drag 
force here roughly equals the radiation pressure. At 
500 km the radiation pressure is "2^ of the neutral drag 
and at 1000 km the neutral drag is '10^ of the solar 
pressure). The radiation dosages acquired during orbit 

( 2 ) 

raising apparently pose no threat to the kapton sail film . 
However the extra shielding needed to protect logic circuits 
would lower the useful payload capacity. 

In the vicinity of earth the sail will encounter 
various plasma and electromagnetic environments. The normal 
precautions that are taken to protect any spacecraft 
during charging events associated with magnetospheric 
substorms will be required for the instrument package. These 
include careful grounding and electrical isolation, and 
cautious use of dielectrics on the spacecraft surface. We 
shall attempt to point out in the following special 
problems associated with the solar sail configuration. 

Most of the time the spacecraft will proceed through 
the very low energy plasmaspheric plasma of density 
•103/cm3 at low altitude to l-IO/cm^ at GSO and above. 

Due to photo-emission from the front surface aluminum and to 
the low energy of the plasma electrons the spacecraft will 
normally be held at some tens of volts positive (positive 
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potentials as high as 100 volts have been observed at very 
high altitude "10^ km). During eclipse the maximum negative 
potential will be kT^ <-36 volts. Neither of these 

conditions are serious for spacecraft operation. During 
substorms this low energy plasma is swept away and replaced 
by energetic plasma clouds whose energy is in the range of 
10 keV with density N'lO/cm^, with associated particle fluxes. 

r = 1 . 7 x 10 ^^/cm^-sec 
e 

r = 4x10®/cm^-sec . 

P 

In the sunlight the photoelectric flux = 2.6x10^°/cm2-sec 
lowers the magnitude of the spacecraft potential substantially. 
With only these three fluxes taken into account and recalling 
that photoemission occurs only on the front side the spacecraft 
will charge to 
<f, = in 

We find in the sun fj’-2kv and in the shade ^^'■-36kV. Back- 
scattering and secondaries reduce these-, typically by a 
factor of two or three. Usually the maximum charging in 
elipse is to approximately kTe, ie -10 kV , and in sunlight 
~-.5kV. Some additional modification of these numbers 
will result if the finite value of the electron gyroradius 
is taken into account. If some dimension of the sail -is 
• small compared to the gyroradius and the plane of the sail 


I 


?! 

m. 
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lies parallel to the magnetic field certain classes of particles 

can 'miss it altogether, thus lowering the electron flux 

and 

ye are assuming that the front and rear of the sail 
are shorted together as discussed in Section 3 . It will 
In addition be necessary to insure that the sail and 
payload are also shorted together. If this were not the 
case they might float up to different and unpredictable 
potentials due to their differing material and geometric 
characteristics, causing problems for electronic systems 
connecting them. 

An important concern for the solar sail is degradation 
of the thin A1 and Cr coatings. The 1-20 kY proton fluxes 
in the magnetosphere are responsible for the sputtering of 
these metallic surfaces. One might imagine that the energy 
increase on attraction to a charged satellite could increase 
the yield. The data we have found however indicates that 
sputtering will pose no special problem. The sputtering 
yield S of protons on silver holds constant at about .05 
atoms/proton over proton energies ranging from 2-20 keV. We 
assume that 10 keV fluxes rp'4x108/cm2 -sec are present over 
some 2 % of the six months the sail spends in earth orbit, 
(corresponding to substorms occurihg in the midnight-dawn 
sector). Assuming a value S'.l atom/proton we find the 
total number of sputtered atoms 

N'1.5x10i“/cm2. 
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The nass of alumintm removed is “7*5x10 ^gm giving 
for the layer thickness removed .3-^* This appears safely 
below the thickness of both the A1 and the Cr . We are 
not certain about effects on the A1 surface reflectivity. 
They are probably minimal because 10^ ^atoms/cm^ represents 
only about one atomic layer. 
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Electromechanical Effects 

We have also considered a number of possible electro- 
mechanical hazards that could conceivably result from 
operations near earth. It appears that none of these pose 
a threat. Most of these effects are associated with the 
passage into or out of eclipse. At GSO it takes some 
four seconds for the heliogyro to pass completely into or 
out of the earth's shadow. We can imagine a pair of 
opposing blades,, one in the sun and one in the shade for a 
period of two minutes. We consider the following phenomena 
that produce mechanical forces that might be detrimental to 
the spacecraft. 

1) Forces due to different charges on the blades; to 
bound possible electrostatic forces we imagine one blade 
in sunlight and one in shade. A voltage of 10kV on a 
blade implies a charge on the blade Q-RV “ (6.25x10^cm)x 
(33 stat volt) = 2x10^esu. If there are two charges of 
this magnitude whose separation is the order of the blade 
length R the force between them is only 10 ^newtons. Thus it 
appears that electrostatic forces on the sail arising from 
spacecraft charging could not possibly compete with the blade 
tension y80 newton. 

2) The current that could flow through the hub section 
during a charging-eclipse event might be on the order of 

I it = = 3x10- %mp. 

This is much smaller than the induction current estimate of 
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Section 7 and thus its effects should be negligible. 

Two types of non-electrostatic forces are apparent. 

These are 3) the torque t arising from lack of sunlight 
pressure on the blade in eclipse and unbalanced 
momentum arising from photoemission from only the sunlit 
blades . 

3) To estimate this effect we imagine a system of 
two opposing blades one illuminated and one not, and assume 
that the moment of inertia of the system is given simply by 
the mass and extent of these two blades. The change in 
angular momentum is aL=t • (A t'2sec) . The torque is provided 
by the radiation pressure P on one blade 

T=Pg J ^rdrw = wPg = 7. 8x10^dyne-cm, 

giving AL = 1 . 5x10^ ®dyne-cm-sec. The moment of i nert ia of 
the two blades is 

C H 

I = 2n r^dr = | nR^ 

where n".4gm/cm is the -linear mass density. We find 
I=6.5x10 ^^gm cm^, and thus a rotation rate 

0 = = 2. 3x10“^ rad/sec = 1.5x10"s degree/sec. 

It appears that this torque would not be noticed. Another 
way to examine the effect is to calculate the Av imparted 
to a single blade and compare to the orbital velocity at 
GSO v~3x105cm/sec. We find av = FAt/m~.2 cm/sec. 

The force due to recoil of photo-electrons emitted 
from a blade at high potential • MOkV) works out to some 

•j.lfxlO "^dyne/cm^j or less than of the solar pressure) and 
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likewise its contribution to the unbalanced force. 

5)' Finally we consider the possibility that the electric 
field associated with some magnetospheric wave disturbance 
could disrupt the spacecraft. The type of wave we are 
interested in will have a period on the order of that of the 
flapping modes of the heliogyro blade, i.e. ~100 sec. Waves 
in this range are known as Pc4 waves and are observed in 
coherent wave trains lasting for some 10-20 cycles. The 
electric field strength can be as high as E".1 volt/meter = 
3.33x10‘^stat volt/cm. 

For a simple model we imagine again a single blade 
charged up to lOkV, so that the charge estimate of 1 ) implies 
a charge density o“4x10~2esu/cm2 . This charge density and 
field imply a force on the blade 

aE~1.2x10‘‘7 dyne/cm3<<Pg, 

and thus the whole spacecraft is not moved significantly. 

We next examine what might happen when this small force is 
applied in resonance with a blade normal mode. We shall 
calculate with a mathematically very simple model. We 
imagine a blade tied at both ends under tension T=385 newton 
3.85x107dyne (or about half the blade tension), the average 
tension along the actual blade. We also neglect damping 
(which is small: Q-1000^^^) and assume E is normal to the 
blade and in perfect resonance with its lowest mode. The 
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wave equation for the blade displacement is 




at 2 


- T 


3X^ 


= F 


ext ’ 


The lowest mode of the unperturbed blade is 

ttX 


u = Sin . Sin wjt 


with u 


= - l~^= c - = 4 
1 RY n ° E 


9x10“2rad/ sec 


RV Ti '' E 

(this gives a period of 127 sec). We also assume for 
mathematical simplicity that the external force has the 
constant shape sin ^ over the blade. In fact the E 
field wave lengths we are considering are much longer than 
the blade and the force is effectively constant in space 
and sinusoidal in time. We are then >considering the 
coupling of this force to the lowest blade mode. The 
wave equation is thus 

fl# = 8 Sin f Sin 

The constant B is given by the force per unit length divided 
by the mass per unit length. 

B = I = iJffi = 2.66x10-fa 

n n sec^ 

'Assuming the displacement is zero at time zero we find 

u(t) = sin ^ (sincoi t -uit cos uit), 

and the amplitude grows linearly, 


A(t) 


-it 
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After twenty cycles or 2,5^0 sec A~7»0 cm. This appears 
safely small. In particular it should be recalled that' to 
compute the force we have used an estimate for the charge 
density that is rather extreme, and also we neglected 
damping. In any case if the heliogyro were developed for 
use within the magnetosphere it would be prudent to 
calculate more realistically the interaction with probable 
magnetospheric wave modes, using a correct heliogyro blade 
model. 

•In siunmary we have found no special problems of a ■ 
plasma or electromechanical nature that would prevent the 
use of the solar sail near earth. 
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5. DEPLOYMENT 'problems DUE TO STATIC ELECTRICITY 


It is important to know whether charge residing on 
solar sail material in the rolled up configuration will 
present an obstacle to its unrolling. The major ’feature of 
the solar sail that seems to minimize such problems is the 
fact that it is coated on both sides with conductors and 
the two coatings are shorted together. Separated static 
charge can reside only on good insulators. It seems likely 
that with careful handling the only surface insulator 
present will be the oxide films that inevitably develop on 
fresh vapor deposited metal surfaces. These are typically 
very thin (30-T00K) and will not sustain a high voltage. 

The resistivity of the film is a few orders of magnitude 
less than that of the bulk material and breakdown potentials 
are of the order of a few volts. It thus appears that we 
will be seeking and attempting to quantify rather small 
effects. This is borne out by the observation that static 
electricity apparently ceases to be a problem after the metal 
film is deposited, In the following we shall consider 

a model of electrostatic forces that should exist even 
between metallized layers, due to the contact potential 
between dissimilar metals, and compare the forces derived to 
that applied in unrolling the sail in the helio-gyro 
configuration. The same consideration should apply to the 
square sail geometry for the case where storage is accomplished 


by rolling. 
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MINIMUM DEPLOYMENT FORCE 

The -un-f-urTi-ng- of th'e' helio-gyro is accomplished by 
spinning the vehicle with thrusters to provide sufficient 
centrifugal force to unroll the blades to a length of 
'150 meters. At this point the blades are given a pitch 
and solar radiation pressure provides the rotation force and 
the blades are slowly released. The final length is 
R “6.250 meter = 6.25 x lO^cm. The rotation angular 
frequency fi during deployment is such that the product 
Rq2 _ constant = Since the rotation period at maximum 

extension is 200 sec. we establish that B^= 6 . 6x10^cm/sec^ . 
The maximum blade tension is specified to be newton. 

We are interested in the minimum force in order to establish 


that it is sufficient to overcome any drag due to static 
electricity. The centrifugal force due to the blade” 
deployed to radius R is 


= "J! 


dr 


r^q^ (R) _ ^2 




where ^ is the 


linear mass density. Thus F scales linearly with R and we 


find F .~20 newton, 
rain 


STATIC FORCE MODEL ^ 

Figure 5.1 illustrates the geometry in the rolled up 
configuration. Opposite charges on the opposing surfaces 
will attract, and retard the unrolling. Free charge could 
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have been deposited on the metallic oxide layers during 
production and s.towage-. We feel that this would have 
dissipated due to leakage conduction by the time of 
deployment. Using values ^ , Table 5.,1 for the bulk AI 2 O 3 
conductivity, which are typically much larger than that 
of the thin film^^^ we find that the characteristic time for 
dissipation of static charge ’ (e~ 10 ), are sufficiently 

small, even at room temperature. 


T(°K) 

p ( flcm) 

o(sec“ ^ ) 

T^(sec) 

287 

Tq16 

9x10‘5 

8. 8x10 3 

573 

3x10^3 

3x10-2 

26.5 

1073 

3-5x10^‘* 

2.6x10^ 

3.1 xIO' 


Table i 1 - Bulk conductivity and charge dissipation time 
for AI 2 O 3 . 

If there is a potential difference between the Cr and A1 
surfaces, however, static charge will remain. Such a 
potential difference might arise due to the contact potential 
between Cr and Al, which is equal to the difference 
between their work functions and the order of 1.0 volt. This 
establishes the amount of charge that remains on the ' 
facing surfaces, and determines the force between layers. 
Contact is established by sail shorting and/or surface 
pressure contact. This effect is well known^^* as a 
mechanism for transferring static charge between dissimilar 
conductors. The geometry is essentially that of a parallel 
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plate capacitor, Figure 5.2. The electric field is E=V^/d, 
surface charge o= E/Mir=V^/4Tid» and force per unit area 
F=oE=V^^/^iTd^ . Appropriate values of d are discussed below. 

The force retarding unrolling will be taken to be the 
average energy per unit length required to separate the 
layers. The metal surfaces come into close contact only over 
a very small fraction of the total area, and, since the film 
is flexible, areas in close contact can be separated by a 
'peeling motion; wheras separating them in the perpendicular 
direction would require much more applied force. We estimate 
■ the forces according to two models. In the first we assume 
that the surface static charge is free to flow as the 
surfaces are separated. In this case the force per unit 
length when the separation is x is F(x) = and the 

average force is 

F„=W / ''c2 dx = w c , 

^ ^d mr 3F 

where w r 800 cm is the width of the blade. We have 
calculated this as the energy necessary to separate capacitor 
plates held at constant potential to infinity. This is 
correct for the unrolling problem as well because the energy 
required is independent of the mode of separation. The 
second model we use assumes that the original amount of 
charge remains on the surfaces as they are separated. In 
this case the force is much larger. Which case applies will 
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depend on the rate of separation, that is on the characteristic 
time for separation compared to the characteristic time 

for charge dissipation. For this second model we have only 
a crude estimate of the force. The model is illustrated in 
Figure 5.3. We simply estimate the energy it takes to take 
the charge from 0 = 0 to above an infinite plane of 

charge, and divide by (Ti/2)Re. We shall argue later that this 
model grossly overestimates the energy and that a much- 
smaller value of is appropriate. The force one ends 

up with is , much larger than F^. 

To complete the force estimates requires estimates of d; 

for the helio-gyro I^30cm. We consider, starting with the 

smallest, the characteristic spacings involved. The VDA 

(vapor deposited aluminum) thickness is d^=10"^cm, with 

uncertainty dg^=10“6cm. We shall disregard the Cr layer. 

The Kapton thickness is dj^=2. 5x10" ‘*cm with uncertainty 

d =5x10“ "205) . The final dimension we wish to consider 
sk 

is what we shall call the "loose wrap", separation 
djj^ =2.5x10“ ^cm^*'\ a rough estimate, of the spacing when the 
sail is rolled up before deployment. However due to the 
presence of battens lying across the sail some 10-305 of the 
area will be pressed together. Thus the force for will 

probably give the most useful comparison with 
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Figure 5.3. Geometry for force at constant charge 



The distances and forces are displayed in Table 5.2. 


d^ (cm) 

F^(nt) 

F^(nt) 

"sv^ 10“^ 

7.07x10-3 

2.12x103 

d^=10-5 

7 . 07x10“** 

- 2.12x103 

dg^rbxIO-S 

1 .i^xlO-** 

85 

d)5=2.5xl0"“ 

2.8x10-3 


'^e«'2-5x10-3 

2.8x10“3 

3. ^x10“2 


Table 5 . 2. Electrostatic Forces Retarding Unrolling 
of Sail., 

We see that if is applicable no problems will arise. 

However, if F^ is involved the forces are larger. We do 

not feel that F_ will apply for the following reasons. 

In the first place raising a charge a distance R implies 

a 6000 volt potential difference. However since the thin 

oxide coatings will only sustain a few volts before breaking 

down(7), they will only separate to a few times d before 

breaking down and thus the force should be of the order of 

a few times F . Aside from this the enhanced conductivity 
c 

for thin oxide films, which is described in the literature 
.as being something like a few orders of magnitude should 

V 

be sufficient to insure that the charge equilibration time 
is short compared to the time of about 60 seconds it takes 
the sail drum to rotate through ir/2 radians. Values of o 
inferred from the literature (eg. fox an oxide 

layer of order 10"^cm thick and .01 in area, (Ref, 8) give 
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)) indicate that this is the case. 


T ^-10* ^sec 

c 

Thus without further detailed informati.on. we tend- to . think that 
the forces are applicable, and static charge should be 
no problem. We should emphasize however that thin film conductivity- 
values vary widely with temperature and method of preparation. 

In summary we feel that there is no strong indication 
of problems pf deployment due to static electricity. However 
we have based our analysis on the average forces involved. 

We have only crudely analyzed the possibility that small 
localized areas that are in very close contact might stick 
together, and the possibility the pliability of the Kapton 
would allow the average separation to decrease. 
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6. ELECTROMECHANICAL OSCILLATIONS OF THE SOLAR SAIL 

The goal addressed in this section is to ansver the 
question "Does the interaction with the solar wind plasma 
cause the sail to develop unstable modes of oscillation?" 

In sailing parlance "Does the sail luff?" This question is,-, 
deceptive in its simplicity hiding a myriad of possible 
mechanisms. To illustrate, consider the equation of motion 

I 

for a membrane 



where for the case at hand 

f = oE^ + P + S 

= electrostatic ^ plasma photon 

pressure pressure pressure 

We are 'interested in perturbations from equilibrium caused 

by the first two terms. As a worst case we shall neglect the 

photon pressure which would appear to be universally stabilizing. 

Our treatment consists of examining a few cases which are amenable 

to calculation then making an informed guess as to the magnitude 

of the problem. This section, while only a few pages long, 

required almost as much research effort as the rest of the 

report. There is almost no directly relevant literature and so' 

much of. the work must be original. There remain numerous loose 

ends and the final answer is far from definitive. Nonetheless, 

with these substantial qualifications, we find that the answer 

to the question at hand is "Probably not." 
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Gas Dynamic Coupling 

- We first illustrate the p.ro-blem by cons-ideri-ng the 
coupling of oscillations on the surface of the solar sail 
to sound waves in the surrounding gas. This ease is amenable 
to calculation and demonstrates one important dimensionless 
parameter associated with electromechanical oscillations. 

The result of the calculation is quite clear. This mechanism 
cannot cause luffing. 

The problem is illustrated by analogy. Consider a 
cylinder, of infinite length, which has a thin membrane (a 
drum head) stretched across it as shown in Figure 6.1. 

Morse and Ingard ^ have considered the similar problem of a 
vibrating string coupled to sound waves in a gas. Review 
of their work indicates the complicated nature of these 
problems. In this section we demonstrate some of the salient 
features of such a system. A simple model is used which 
substantially reduces the algebraic complexity of the problem. 
The solutions presented here are more easily imderstood than 
those in Morse and Ingard and exhibit substantially the 
same physics. 

If the drum head is struck with a hammer it will 
oscillate with a characteristic frequency determined by the 
parameters of the system. The oscillation of the drum head 
will set up sound waves in the cylinder. These sound waves 
will propagate down the cylinder away from the drumhead. 
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membrane-^ 

V = (T/p)^ >>Cs 


Figure 6,1. Waves propagating on a membrane coupled to a 
surrounding gas. When the wave speed on the 
membrane V = (T/p)^ << C g» the sound speed in 
the gas, then the mass of the gas loads the 
string. When V >> Cs sound waves are radiated 
from the string in a process analogous to 
Cherenkov radiation. 
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Consider for the moment the modes of oscillation with dimensions 
small. -compared- to that oT the cylinder i.e. wavelengths X<<R. 

In this limit one may neglect the details of the boimdary ^ 
conditions and treat the membrane as an infinite sheet embedded 
in an unbotinded medium. What are the modes of oscillation of 
this system? 

Waves propagate on an infinite membrane with the 
characteristic dispersion relation 

where w = 2iif is the characteristic frequency 

K = 2ti/x is the wave number for X the wavelength 
T is the tension in the membrane 
P is the density of the membrane. 

The characteristic dispersion relation for sound waves 
in a gas is 

0) = C^K 

where C Is the sound speed which depends upon the 
s 

temperature and molecular composition of the gas. 

This is not a complete description of the problem at hand. 
The presence of the gas around the membrane causes additional 
forces to act upon the membrane thus modifying its characteristic 
dispersion. Similarly the dispersion relationship for waves 
propagating in the gas is modified by forces exerted by the 
membrane upon the gas. The oscillations of the membrane and 
the gas are thus coupled. 
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Two limits suggest themselves. Whenever the wave speed 
on the membrane is less than the sound speed in the gas, i.e. 
(T/p)^<<C the gas can respond essentially instantaneously to 

3 ^ 

deformations of the membrane. Thus an oscillation propagating 
slowly across the membrane will simply push the gas back and 
forth. In the alternative limit (T/p)^>>Cg the gas cannot 
respond instantaneously to the oscillations of the membrane. 

In this limit an oscillation propagating across the membrane 
will excite waves in the gas which propagate away from the 
membrane. The two regimes are illustrated in Figure 6.1. 

- These effects will now be illustrated by detailed calculation. 

The membrane is described by the linearized wave equation 





where n is the displacement of the membrane 


in the y direction 

and p+ are the pressures exerted by the gas 
for y< 0. 

The gas is described 'by equations for the conservation of 
number density n and momentum, and an equation of state 
relating the pressure p+ to the number density. These are 
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p = 




n 


n 


where m is the mean molecular weight 

Y = 1 for isothermal expansion 

Y = 3 for adiabatic expansion 
Combining the gas equations one finds 


. C 2 

3t^ S 


SX 


9i4 

9y^ 


= 0 


which describes the propagation of sound waves in the gas 

, - ... M, 

s 


characterised by a sound speed C_ = (yP /n m)' 

0 0 


(yT /mV 
0 


where T is the temperature. 

The gas and membrane are further coupled by the requirement 
that the component of gas velocity in the y direction match 
the velocity of the 7 displacement of the membrane. 

Thus ?y = *1^ on the membrane. Using the momentum 
equation we relate this to the pressure and find 


3t^ ' ^0® 

which for linear perturbations of the membrane is to be 
evaluated on equilibrium surface y = 0. 

In the limit (T/p)^<C the waves do not propagate into 
the gas and the solutions for the pressure may be represented 


pe““y 

y>0 

pe“y 

y<0 


where a is a constant to be determined from the gas propagation 
equation 
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Combining these equations one finds 


p + = + 



p.-p+ = 



Thus we find that the equation describing the coupled 


oscillations is 


P + 


2n m 

O 




This equation admits to a very simple interpretation. 

The membrane mass density per unit area p is replaced by an 
effective mass density per unit area p* = p + (2n^m/a)* So 
we see that the mass of gas contained in the perturbed pressure 
sheath of thickness L = 1/a simply loads the membrane. 

Now let us examine this limit in slightly more detail by- 
considering a specific perturbation of the form 
=. n sin(kx-ti)t) 

/ pe’“^sin(kx-(i) t) y>0 


P = 


pe°^sin(kx-ut) 


y<0 


The dispersion relation becomes 

l^p ■+ 0)2 - T K? = 0 

where — 

•tj ? ^ ^ - 


This dispersion relation is cubic in thus exhibiting three 
double modes corresponding to the modified sound waves on either 
side of the membrane? and the basic membrane mode. 
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Now let us examine the effect of a gas flowing across the 
membrane. While a number of additional complications could 
be introduced, we shall only examine the effect of doppler shift 
which causes the frequency associated with the moving medium 
to be transformed from w tou • Thus pressure perturbations 

associated with a gas moving with velocity v = v e^ will obey 

O 'Ww 

the equation 


70 



(o)-K V )2 - K2 c_2 
0 


p + 1^ = 0 


which has evanescent solutions 


P = 


/Pp“oypi(Kx-u)t) 




y->o 


y<o 


whenever 

a2 = K2 


1 - 


KC. 



>0 


• Now consider the case ^‘^<V <C_. The appropriate 
solutions are the evanescent ones since a2>0. The boundary 
condition is unmodified but the momentum equation 
is doppler shifted. One obtains 
1 / 1 \ _J_ !£■ 


ri - 


u I 0) 


-KV 


1 

n m 
0 


ay 


and hence 


p-"p+ = u(w-KV ) 


= - t]j2 


(i¥j. 



So there results the dispersion relation 

m\‘ 

-^j co2-T K2 = 0 

effective mass p* = p- (KV^/w) (2n^m/o^ may 
become negative in which case an instability would develop. 
The physical mechanism is as follows. The mass displaced by 
the sail is convected along by the flow thus shifting by 180 
degrees the phase with which the mass loading is applied to 
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the membrane. The resulting instability grows by extracting 

free energy from the flow. 

A basic dimensionless parameter e = (n^mL/p), where 
L = is a typical length, characterizes the gas dynamic 
coupling. This parameter simply measures the effective mass 
loading which the gas applies to the membrane. 

Tot the solar wind n^m < 1 .6x10”^^gm/m^. For the sail 
P = 3 gm/m^ thus e = 1 requires L = (p/n^m) > 2x10^^ meters. 

This astronomically large number indicates that gas dynamic 

I 

mass coupling is not significant, since L is many orders of 
magnitude greater than a typical sail dimension. Alternatively 
note that in the streaming case L = (V^/tu) where = 800 Km/sec 
and w = 5x10” 2 rad/sec are typical for the solar wind and the 
sail respectively. Thus L = 1.6x10'^ meters and hence e = 10 
so we are many orders of magnitude away from insta bili ty. 

Thus we can say with certainty that the solar wind is of 
such low mass density that it cannot couple significantly to 
the sail through these mechanisms. 
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Resonant Forcing 


Obviously we must look elsewhere for potential sources 
of instability. Treating the solar wind as a simple gas 
rather ' than as a plasma neglects the important electrostatic 
couplings between the solar sail and solar wind. To better 
understand these couplings we now examine the case of resonant 
electrostatic forcing of the sail. These are forced oscillations 
caused by those fluctuations of the solar wind which are of 
such a frequency as to resonate with the natural frequency of 
oscillation of the sail. 

Consider a membrane subject to external forces. The 
equation describing linear oscillations of such a membrane is: 



where AF is the perturbing force exerted on the surface of 
the membrane. In the previous section we examined perturbing 
forces caused by a pressure differential across the membrane. 
Now we examine electrostatic perturbing forces given by 


aF = o^AEy + AoE^q 
E 2 


where U = (E ^/l+n) is the electrostatic energy density 

^ y 

evaluated at the membrane surface. 

In equilibrium, as has been discussed in earlier parts 
of this report, a plasma sheath approximately one Debye , length 
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thick will form about the sail, A typical particle will transit 
this sheath in times t which for ions is typically a 

few microseconds. The natural frequencies of oscillation of 
the solar sail are many orders of magnitude less than this 
thus the sheath responds quasi-statically to fluctuations 
in the solar wind. Therefore we take 

AU = U ( Aa/ a ) 
c c 

where U^'v^nm^U^^'vIO ^erg/cm^j which is the wind streaming 
energy density, and (Ao/o) is the relative amplitude of a 
fluctuation in the solar wind. 

Typical \anforced oscillations of the membrane are of the 

form 


fl = n sin(ux/L)sina>it 
where = (tt/L)(T/p)^ 

Now if we consider a resonant fluctuation for which 

— (x,t) = — sin(7Tx/L) sinajjt 

then the sail will respond with a secularly growing oscillation. 
For an initial displacement of zero, one finds the solution 


T1 - 




I —^1 ( 

7^ V p / i ® / 


And thus the amplitude 


sin(Trx/L) (sinioit - a>it cos uit) 
of the secular growing term is 



Now numbers 

= 10“‘^ergs/cm3 


10"8 joules/m3 
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p = SxIO’^Kg/m^ 

u)i= 4.9x10‘2 rad/sec 

which result in an amplitude 



Thus (t^a/a) 0.5 5 indicating a fifty percent change in 

amplitude of the solar wind, would require uit 33? i.e. 

33 cycles, in order to produce a wave with amplitude A 'v- 1cm 
on the surface of the sail. Such long trains of large and 
coherent fluctuations in the solar wind are quite unlikely. 
Thus we conclude that resonant forcing should not be a 
serious problem. 

To gain insight into this r.esult rewrite the expression 


for the amplitude as 


n'v 


^ I II Aa I /til 1 tV 
n V T ylct /(2 ti ) 


Noting that in equilibrium 


-T 




where U , is the photon energy density. Hence by dimensional 


analysis one has the gross estimate: 

T-v 


n-. 


And therefore 




The significant term is (U^/U^)< 10"^ and thus typically many 
cycles are required before a large oscillation can build up. 
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So once again we see that the sail does not luff. The 
dimensionless parameter this time involves the energy density 
ratio, CU L-^/Tn C'U' /U ) which is quite small. 

' E o e V 

Debye Sheath Coupling 

An alternative approach to the electrostatic perturbing 
forces is next examined. Again we consider a membrane perturbed 
by an external electrostatic force, and hence governed by the 
wave equation. 



This time we take isF to be caused by self-consistent 
perturbations of the Debye sheath. The rippling of the sail 
perturbs the sheath, which in turn causes perturbing forces 
to act back on the sail. Hence this calculation is close in 
spirit to the earlier calculation of gas dynamic coupling. 
Now, however, we are examining coupling to the plasma Debye 
sheath. 

We take a simple exponential model for the equilibrium 
sheath. Forces from either side of the membrane add 
vectorially thus we present the force derivation for a 
single side and only indicate the result for two sides. 

For the front side, y>0 one has an equilibrium sheath 
0 = 0^exp(-K^y) electrostatic potential 

Eyo = exp^jjY) electric field 

~ surface charge density 

where = (47rn^e2/T)^ 
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The perturbing force is 
AF = AaEy^ + o^AEy 

where AE = 4-nAa are the perturbed electric field and 

y 

charge density evaluated at the membrane surface. Note that 

© © 

3E 

AE = ■ — ^ 6y + 6E , 6y = n 

y 9y y 

where term © represents the contribution to the total 

t 

perturbed electric field which results from evaluating the 
equilibrium electric field on the perturbed boimdary and 
term © represents the perturbed fields evaluated at the 
equilibrium boundary. 

To determine the pertirrbed fields we need a prescription 
for the plasma. Again, we note that for times t long 
compared to an ion transit time through the Debye, sheath 
^ ^sec j, that the quasistatic approximation is 
valid. In this approximation the potential obeys the Debye 
equation 

v2 60-Kj^2 50 _ Q 

which, for perturbations varying as sin kx on the membrane 
has the solution 

2 ^ y 

60 = 60e“ sin Kx 


so that AEy. = 4TrAa 

960 ~^yo 
'9y 3y 


9E 


6y 
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We examine the case of a conducting membrane, which is 
an ap.propri-a-te- representation of the sail. (For wavelengths >> the 
sail, thickness the two sides of the sail move together and no 
charging of the dielectric takes place in first order.) Thus 
ve demand that the total potential, which again consists of 
two terms like (3) and @ be zero on the membrane. So one 
has 

0 = 6y + 60 = 0 

Hence using the equilibrium results one finds that 

60 = 


and so 

iEy = I I 

The force on one side is thus 

] 'y 

The total force is obtained by simming the forces acting on 
the front and back sides. Thus 




AF = AF 

f,b 


where f = front, b = back 




/ 




6y 


f ,b y2Ti 

and therefore the wave equation may be written as 


P = -TK26y 




Sy 
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in the long wavelength limit K<<K^ and this reduces to 


= . ik2 


8t 




5y 


A plasma which is the same on the front and back will 
exhibit unstable modes whenever 


2 . 


> K^T 


for the solar wind 


U = 
c 


4ti 

'v 10'"^ergs/cm2 


while 


'o'l 0“ ** ergs/cm ^ 

Thus no instability will occur. Again we find a dimensionless 

energy ratio which is quite small. 

Note again that T'v.(L2U /n ) where n is the equilibrium 
displacement of the sail. Thus the stability parameter can 
be written as 

Ka'g’o' 

■ “t’KaT - 

Each element is this parameter is much less than one thus ^ 
is quite small corresponding to a very stable situation. 
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Application to Heliogyro 

The detailed calculations have used a membrane equation 
as description of the basic sail. This is adequate for 
illustrative purposes and is perhaps a reasonable representation 
of "the square sail. A detailed analysis of the heliogyro could 
be based upon the uncoupled small-motion equations for twist, 
vertical deflection and inplane deflection. In view of the 
strong indication that luffing is not a serious problem, which 
has been derived using the membrane equation, the use of the 
more elaborate equations seems unjustified. In the event 
that such analysis were necessary one would again add the 
electrostatic and gas dynamic forces to the basic equations. 
Concluding CoTninents 

These calculations indicate two dimensionless parameters, 

= (n^m/pa) the ratio of mass densities and 
the ratio of plasma energy density to photon energy density. 
Both are much less than unity and thus the gas or plasma 
forces acting upon the sail are much less than the forces 
caused by photon pressure. No mechanism has been found which 
would significantly amplify the impact of the plasma forces 
thus no instabilities develop. 
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7. INDUCTION CURRENT AND INDUCTION FORCE 

When a conducting satellite moves thrgugh_^a stationary 

magnetized plasma an electric field E = - ~ is set up 

within it, due to the polarization of the free charge 

produced by the motion perpendicular to B ’ , 'Figure 7.1. 

Likewise since B field lines can be regarded as frozen into 

the solar wind plasma motion of the solar wind plasma pas.t 

^ _ y X B 

the satellite with velocity V. produces the field E - _§ . 

s c 

The potential variation on the satellite surface causes 
electrons to land preferentially on the more positive 
portions, and protons to land preferentially on the more 
negative portions. Also photo-electrons emitted from the 
positive parts are more likely to be attracted back to 
the surface. 

1 

This situation causes a steady current to flow through 
the satellite. The current is calculated by using the 
potential at each point on the surface to calculate the 
plasma current into the satellite. The potential is a linear 
function of the dimension x in the Vg x B direction. The 
requirement that the net current into the satellite equal 
the net current out of the satellite gives the potential. 

The possible consequences of this current in our view are 
twofold, the current crossed with the magnetic field produces 
forces and torques , and the current dissipates energy in the 
aluminum sail film. We shall argue that the induction forces 
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are far smaller than the radiation pressure, and that the 
energy dissipated by conduction is negligible in comparison 
to the sunlight energy absorbed by the aluminum. 

It is not necessary to perform a detailed calculation 
to show that the above effects are small. Figure 7.2 shows 
the general situation with the currents. For this diagram 
we assumed that the potential changes sign somewhere along 
the blade. -This is not necessarily true. The net current 
per unit area into or out of the sail at any point x is' 
given by the sum r(x) of the three fluxes at that point. 

A crude bound on the maximum current along the sail can 
be made by taking the largest of the three current densities 

= r , = 2.6 X 10ii/cm2- ser (at .3 a.u.). multiplying 

e ph 

by the electronic charge and one half the area and assuming 
this current flows in one half the blade and out the other . 
half. This current is 

I = . • (w=8x102cm)(R=6.2J4x105cm) = 3* T2x10^°esu/sec, 

max ph 

I - = 10.5 amp . 

max ^ 

This total current multiplied by the length and 
magnetic field and divided by the blade area gives a value 

for the total induction force per unit area 

pi *R • (B=58x10"Sgauss) 

F = max j 

c*R‘W 

F s 7.54 x 10~'^dyne/cm2 . 

This is about .15? of the solar radiation pressure at .3 a.u. 
and thus seemingly negligible. 
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(a) Potential on Blade 





Figure 7.2. Sail potential and particle fluxes. 
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The current density in the VDA implied by is 

j = ^max = 3.9 X lO^^esu/cm^-sec . 

w(t=1000S) 

Using the approximate value for the conductivity at 525 K 
o~2 . 6x10^ "^sec~ ^ we find the volume energy deposition rate 

^-£ 1-5 B5x10"^^^^^^f = 5.85 watt/cm^. 

dt'o cm-^-sec 

This implies a per unit area energy absorption of 5.85 x 10 ^ 
watt/cm^. The VDA absorbs about 10^ of the incident optical 
energy, or at 0.3 a.u. about 0. l4 watt/cm^. ^ Thus the 
induction current should not affect the thermal balance of 
the sail. 

In Section 9 we point out that the induction current 
might possibly interfere with magnetic field measurements. 
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8. EFFECTS OF SOLAR SAIL OPERATION ON THE PAYLOAD 

In this section ..ve consider the effects of the presence of 
the solar sail on the operation of the payload itself. . 

The possible areas of concern are 

1 ) Safety 

2) Effects on particle measurements, 

3) Effects on field measurements 

4) Electromagnetic Interference (EMI). 

We will not consider the possible effects. on various optical 
observations since this will depend on the specific configuration. 
We do note however that the periodic illumination of various 
parts of the payload by the rotating blades could cause spurious 
readings. Likewise light scattered from matter emitted from 
the sail could influence readings of very low light intensity 
measurements . 

In the proposed configuration, the payload is asymetric 
about the sail spin axis. Further, one part of the payload 
is shielded from the sun with a large shade. Several aspects 
of this proposed configuration could have serious effects on 
the operation of the' scientific instruments as well as perhaps 
presenting a hazardous situation. 
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Safety . Because the payload part of the solar sail vehicle 
is probably composed of materials with different surface 
characteristics from the sail itself, the equilibrium electro- 
static potential that it would assume in the undisturbed 
medium could differ considerably from that ass.umed by the sail. 

If the ground between them depended on a potentially intermittant 
connection through the despin mechanism, then discharges could 
occur with large currents being injected Into the payload proper. 
Therefore we emphasize the need for making firm connections 
between spacecraft chassis ground and the sail. The current 
through this joint should be of the order of a few tens of 
microamperes and be more or less constant. 

The alternative of completely isolating the payload seems 
to be less desirable for several reasons. It would make 
communication between the parts more difficult (i.e. optical 
couplers needed), and the varying differential potentials 
would affect adversely particle measurements while the payload 
is still attached to the sail. 
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Effects on Particle Measurements . The presence of the solar 
sail can influence the attempts to measure the particles in 
the undisturbed plasma in several ways: 

1) Production of extraneous particles 

2) Creation of local electric fields, and 

3) Creation of local magnetic fields. - 

As has been discussed in other sections of this report, the 
sail will produce copious amounts of photoelectrons, back- 
scattered electrons, and secondary electrons. Smaller amounts 
of neutral particles and positive ions are also produced. In 
the most likely charge state, the total spacecraft will also 
accelerate (by the potential of the whole structure) either 
the electrons or ions while retarding the other species. Since 
there is most likely no way to eliminate all of these effects , 
the experimenters will be forced to design instruments which 
operate in the sail environment. For instance, an instrument 
designed to measure low energy (less than 20 eY) electrons in 
the vicinity of a comet will risk excessive counts on its 
sensors if operated for long periods in the cruise mode. There 
fore the instrument must have command modes available to avoid 
this energy range during cruise, or only be operated for short 
test periods. 

As has been shown by "Whipple^^^ using data from ATS-6, 
careful analysis of particle data as a function of energy and 
angle with respect to the spacecraft coordinates can allow the 


88 



experimenter to deconvolve the effects of local electric 

fields and thus sort out many of the locally produced 

{ 

particles from the natural ones. However this analysis is 
time consuming and costly. 

No algorithm for this kind of study has yet been 
developed for the solar sail, but the nature of the heliogyro 
might simplify the work somewhat because of the periodic 
modulations to be expected as the gyro rotates. 

Creation of local electric fields by the sail will distort 
the trajectories of the incoming ambient particles. These 
fields are also responsible for reflecting the locally 
produced particles back on paths which allow them to be counted 
by the payload instruments. ‘ As was alluded to above, some 
mapping of these fields can be accomplished by studying the 
measured distribution function of all particles reaching a 
given detector. This might be the only way to measure such 
fields if long booms cannot be deployed from the payload 
before jettisoning the sail. (We do not consider here the 
interesting posibility that the gyro blades themselves can be 
used as electric field probes.) 

These locally produced fields will most likely consist 
of an overall charge, a potential minimum produced in the 
plasma sheath, and localized irregular fields due to differential 
charging. But to this list we must add another possible 
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candidate due to the proposed configuration. The various 
opti.c.al sensors are- sh-i-elded from the sun by a large shade. 

This will cause dielectrics in the area (such as the lenses 
themselves) to charge more negatively than the main body. 
However light reflecting from the blades of the gyro as they 
rotate under this section can produce enough photoelectrons 
to discharge the dielectrics periodically. This could 
produce a large modulation on local electric fields which would 
be very difficult to untangle in the particle data. Similar 
comments can be made about any shadowed portion of the 
payload, but this effect is surely heightened by the presence 
of the sun shade. 

As disruptive as these effects will probably be on 
particle measurements, the effects of locally produced 
magnetic fields might be worse. If we assume an induction 
current in the blade of the order of 10 amperes (see Section 
then we can easily estimate the magnetic field to be expected 
at the instrument aperttire due to this contaminant. For 
probable dimensions, the induced magnetic field i's in the order 
of hundreds of gammas (lO'^gauss = 1 gamma). In the solar 
wind the undisturbed field in typically of the order of gammas 
or less. 

The effect of this field will not be noticable on most 

ions or on electrons of more than a few hundred electron volts 

energy. However, the lower energy electrons will have their 

trajectories so distorted that one might not be able to 

deconvolve the bending the underwent in this magnetic field. 
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Naturally, this problem can be solved much more easily 
than the problem of the stray electric fields. All one has 
to do is to prevent the existence of such large currents near 
the payload. This should be possible without compromising 
other mission objectives. 

Effects on Field Measurements . Of necessity, local field 
contamination has been discussed in the previous section 
on particle measurements. All that needs to be said here is 
that deployment of long electric antennae is not advisable in 
the presence of the solar sail. Even stub booms might be 
saturated a large fraction of the time. 

The example of the possible levels of magnetic contamination 
given in the "previous section shows that a magnetometer on 
the payload would have a difficult time indeed trying to measure 
natural fields in the presence of a periodically varying 
field which is potentially orders of magnitude greater than 
the desired signal. 

ElectromaEnetic Interference . EMI induced in the payload by 
the solar sail is a possible hazard in addition to the obvious 
safety aspects. Excessive noise could trigger counters, upset 
logic, etc. However, at this time we have been unable to 
isolate' any aspect of the solar sail itself which would 
produce excessive EMI, providing that the grounding and sail 
shorting is carefully done. 
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9. EXPERIMENTAL PROGRAM 

Even th.Qugh mos.t of the anal-y-si-s done in this report 
indicate that possible sail-plasma interactions are an order 
of magnitude or more smaller than other interactions, an 
experimental program should be undertaken to verify the various 

results . 

It would be beyond the scope of this study to prepare a 
complete experimental program here, but a similar program that 
might be of use to the study of the solar sail has been 
prepared by MAYA for NASA Lewis Research Center under a 
subcontract to Systems, Science, and Software (S^) of La Jolla. ^ 
The purpose of the report was to develop plans for the experiments 
which must be conducted in order to verify the ground test 
mathematical model (GTMM) of spacecraft charging which has 
been developed by S^ for Lewis Research Center. 

Special problems that involve testing of the sail involve 
its size. For instance, electrical-mechanical oscillations 
might not scale linearly, and testing a full-sized blade is 
not feasible. 

Therefore we suggest that the best policy is to test as 
large samples as possible by placing them in high-vacuum 
systems with appropriate irradiation by particles and photons. 
VJhile close watch should be kept for induced mechanical 
effects, the prime observations should be the leakage current 
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and induced potentials. Simultaneously EMI produced from 
sporadic breakdowns should be monitored to give an idea of 
what might be expected in orbit. Because of the very thin 

material used, large potentials cannot be generated before 

( 

breaking down. Therefore we expect much less spectacular 
displays with solar sail material than has been seen in 
similar tests at Lewis on normal spacecraft materials. 

Contamination of the environment probably cannot be 
estimated accurately from ground tests since the interactions 
with the vacuum system itself will also contribute extraneous 
particles . 

Some attempt should be made to assess the lifetime of 
the sail in simulated solar wind conditions. We have not 
attempted to assess the impact of induced sputtering or the 
effects of cross-linking of the polymers by radiation, but 
these and other effects which are indirectly related to the 
natural environment could combine to shorten the expected 
lifetime on orbit. 

Finally, we strongly recommend that anattempt be made 
to simulate the deployment of a sail blade under as realistic 
conditions as possible. As can be seen from Section 5 on 
deployment, there is much more than first meets the eye about 
possible problems in unrolling this very thin membrane. Even 
when the front and back surfaces are shorted together, there 
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might be some net force induced by the environment. As in 
the previous case, scaling here is difficult.. -The .expected- 
forces are quite small, but might be significant. A possible 
observation would be discharges at the separation line. 
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